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Abstract 

We present two descriptions of the the local scaling and shape of ideal 
rings, primarily featuring subsegments. Our focus will be the squared radius 
of gyration of subsegments and the squared internal end to end distance, 
defined to be the average squared distance between vertices k edges apart. 
We calculate the exact averages of these values over the space of all such 
ideal rings, not just a calculation of the order of these averages, and compare 
these to the equivalent values in open chains. This comparison will show 
that the structure of ideal rings is similar to that of ideal chains for only 
exceedingly short lengths. 

These results will be corroborated by numerical experiments. They will 
be used to analyze the convergence of our generation method and the effect 
of knotting on these characteristics of shape. 
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1 Introduction 



Long strings of connected molecules, called polymers, are central structures in 
life and physical sciences, as well as engineering. Prominent examples are DNA, 
proteins, polystyrene, and silicone. With regard to DNA, Fiers and Sinsheimer 
first showed that the DNA of a specific virus is a single-stranded ring [8] . Because 
of this closure condition, DNA, like other closed polymers, can be knotted. In 
1976, Liu et al. discovered examples of knotted DNA, which was followed by the 
discovery of topoisomerases, enzymes which knot and unknot DNA [2j [13] . These 
discoveries suggest that knotting plays an important role in the behavior and shape 
of DNA and other polymers. 

We will examine ideal rings, embedded equilateral polygons in M 3 , which pro- 
vide a model for polymers under the ^-condition, where excluded volume can be 
ignored [31 |9l [22] . We will compare ideal rings with ideal chains, random walks 
in R 3 , to determine the effect of the closure constraint as well as the effect of 
knotting. 

These models will be used to analyze the relationship between shape, scale, 
knotting and two physical characteristics specific to their local structure. We define 
the average squared radius of gyration of subsegments, calculated by averaging 
the standard squared radius of gyration of each subsegment of length k, and the 
squared internal end to end distance, the average distance between vertices k edges 
apart. Rather than finding approximations of these averages, we will determine 
the exact theoretical averages for these descriptions of shape, as well as discuss 
numerical simulations and show how these characteristics are affected by knotting. 

In numerical studies examining squared end to end distance, squared end to 
end distance is refered to internal end to end distance or two point correlation J2TJ 
[2H EE] ■ These characteristics are simplified in open chains because subsegments 
of length k in chains of length 100 are identical to subsegments of length k in 
chains of length 100, 000: the ambient length has no effect on the behavior of 
subsegments. However, in the case of ideal rings, subsegments of length k in a ring 
of length n do not behave the same as a subsegment of length A; in a ring of length 
m ^ n. Thus, subsegment behavior relies on both k and n. It will be shown that 
the structure of ideal rings is similar to that of ideal chains for only exceedingly 
short lengths, which may correspond with the difference in ^-temperature between 
open chain polymers and ring polymers [22]. Understanding in the case of closed 
chains may also allow us to use these charactereistics to identify and describe the 
knotted portions of open chains [HI EI] ■ 
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2 An Introduction to Ideal Rings and Their No- 
tation 



Definition 1. An ideal ring, P, is an n-edged equilateral polygon embedded in IR 3 , 
with one vertex at the origin. Let e±, ...,e n be unit vectors such that the kft 1 vertex 

k 

of P is Vk = ^^ e i- Let each e, be called an edge vector of P. The fact that P is 

i=l 

n 

a polygon is equivalent to requiring closure, that is, that v n = ^^e^ = 0. Let V n 

i=i 

denote the space of all such polygons. 

The careful reader may note that these ideal rings are more specific than usual. 
We require these polygons to be based at the origin, and be oriented. This defini- 
tion makes notation concise, and does not affect the averages we will calculate, as 
they are independent of the base point and orientation. 

For each P EV n , the edge vectors are all identically distributed. Therefore, 
when i ^ j, (e« • ej) is a random variable that does not depend on i and j. 

We will compare these with ideal chains, a similar population but without the 
closure constraint. 

Definition 2. An open chain, W , also sometimes called an ideal chain or ideal 
open chain, is an n-edged random walk in IR 3 , where each edge has unit length. Let 

k 

ei, ...,e n be unit vectors such that the kft 1 vertex is given by Vt = ^^ e i- As with 

i=i 

ideal rings, let e; be called an edge vector. Let W n denote the space of all such 
open chains. 

Again, the careful reader will see that these are also based at the origin for 
ease of notation. 

3 Rigorous Calculation of Theoretical Averages 
for Ideal Rings and Ideal Chains 

3.1 Squared End to End Distance in Ideal Chains 

For comparison, we will first calculate the squared end to end distance in ideal 
chains, as in [S]. 

Lemma 1. In W n , the average squared end to end distance is n. 
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Proof. Let W G W n . Calculating the end to end distance of W: 



i=l i=l 

n / n— 1 n. 

= ^( ei -e 4 ) + 2 ^^(e 4 -e,) 

i=l \i=l j=l 

For random walks G W n , the direction of each edge vector is completely 
uncorrelated with the direction of the previous edge vector, and each edge vector ej 
is a uniformly distributed, unit length, random vector. Thus, for i = j, (e^-e^) = 1, 
and for i ^ j the average value of (ej • ej), (e* • ej), is 0. Thus we have that the 
average end to end distance, taken over all W G W n , is 




□ 

3.2 Average Edge Product in Rings 

Now we will consider ideal rings. As noted before, (e^ • ej) is independent of % and 
j. Here, we find an average value for (ej • ej), taken over all i, j for all P G TV 

Definition 3. Consider the space of ideal rings, V n , for some n. Let r n denote 
the average of the set lZ n , 

lZ n = {(ej • ej) : ej, ej are edge vectors of some P G V n , % ^ j}- 

We will call r n the average edge product. 

The following Lemma is generally known, and is foundational to the following 
proofs [TO] . 

Lemma 2. For all n G N, the average edge product, over all ideal rings of length 

n, is given by r n = . 

n — 1 
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Proof. For all P £ V n , = 0. Then squaring both sides we have 

i=i 

n n 

(0-0) = = ($»•($» 



i=l i=l 



= 5Z(e;-ei) + 2 (ei-ej). 

l<i<j'<n 



i=l 



Taking the average over all P £ V, we replace (e; • e^) with (e* • ef) = 



n n(n — 1) 
= n + 2^ J -r„ 



Solving for r n we have 



— n 



-1 



r„ = 



n n(n — 1) n — 1 



□ 



3.3 Squared End to End Distance 

In open chains, the average squared end to end distance of a subsegments of length 
A; in a chain of length n is identical to the average squared end to end distance 
of chains of length k. That is, the ambient chain has no effect on the shape of 
the subsegment. The same is not true in ideal rings, where the closure constraint 
plays a pivotal role in local scale, as we will see. 



Definition 4. For any P £ V n , define 



d P (k,j) 



j+k 

i=j 



j+k j+k 
i=3 i=3 



the squared distance between the andj+k^ vertex. We call dp(k, j) the squared 
end-to-end distance of length k at j of P. Define A(k,n) to be the average value 
of 

{d P (k,j) : 1 < j <n,PeV n }, 

that is, let A(k, n) denote the average squared end to end distance of a subsegment 
of length k in an ideal ring of length n. We will call A(k,n) the average end to 
end distance of length k. 
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Figure 1: In the above image, we have end to end distances marked for k = 2, 
k = 3 and k — 4. We would take the average of their squares to find the average 
squared end to end distance for k = 2,3 and 4. 



Theorem 3. A(k, 



n 



kin — k) 
n — 1 



Proof. First we have that 



dp{k,j) 



j+k 

i=j 



n — 1 



j+k j+k 

i=j i=j 
j+k 

= }Xej • e») + 2 2J ( e * ' e »0 
= + 2 2J (e* • e m ) 

j'<j<m<j'+fc 

We have shown that independent of P, k, i and j, the average value of (ej-e m ) = 
1 

. Therefore the average value of dp(k,j) is, independent of P, i and j, 



k + 2 r n . This average is precisely A(k, n). We conclude that 

j<i<m<j+k 



A(k,n) = k + 2 V" r„ 

j'<i<m<j+fc 

,(*-!)(*) -1 



A; + 2- 
n — 1 



n — 1 
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Therefore the average value of II Vk II 2 is A(k,n) = . □ 

n — 1 

We can verify some key features of this identity immediately. At k = 1 and 
k = n — 1, we have that the average squared end to end distance is 1. Vertices one 
edge apart must be distance 1 apart, as we would expect, so the theoretical aver- 
age agrees with the physical reality. Also, the function is symmetric about |, as a 
segment of length k and its complement, a segment of length n — k, should have the 
same end to end distance. Lastly, the function achieves its maximum at |, as after 
that point, on average, the end to end distance must decrease due to the closure 
condition. Witz showed through a different argument that squared end to end dis- 
tance scales like k ( n ~ k > ; which is close to our value [25]. This is an approximation 
though, as for k — 1, we have ^ 1, though it is close for large values of k and n. 

As in [6] and LemmaJTJ for open chains, the average squared end to end distance 

for ideal chains is n (or nb 2 where b is the length of the segments.) The squared 

end to end distance of subsegments of length k within an open chain of length n 

will be identical to the squared end to end distance of an open chain of length 

k. Therefore, in the open case, the average squared end to end distance of a 

subsegment of length k within a chain of length n is k. As we can see in Figure 

[2j the average squared end to end distance of subsegments of an open chain and 

the average squared end to end distance subsegments of an ideal ring look very 

different for subsegments of very short lengths when n = 50. We observe that for 

the average squared end to end distance of segments of length k in open chains to 
th 

be within ^ of the average squared end to end distance of segments of length 
k in ideal rings, we would need k ^ or smaller. For length 50, as in the figure 
below, k must be length 1, which is a very strong restraint. 

3.4 Average Center of Mass 

For the following definition and Lemma, it is helpful to recall that these ideal rings 

n 

are based at the origin, = 0, as we stipulated originally. 

i=l 

Definition 5. Let n6l Let P € V n . Define Cp, the center of mass of P, to be 



c P 

n 

k=l 



1 

- y Vk, the average of the vertices, 
n 

fe=i 

Lemma 4. For any P £ P n , 

^(f>-A; + l) 2 + 2 Yl (n-J + l)(n-t + l)(e t - ej )). 



cp \\ 2 - 

k=i i<*<i<™ 
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10 20 30 40 50 



Length of Subsegment k 



Figure 2: Here we compare the theoretical average squared end to end distance of 
subsegments of length k for open chains (blue) and ideal rings (red and dashed), 
both with total length n = 50. 



Proof. 

n L — 4 n z — ' 

k=l k=l 
1 n k n k 

^EE^EE^)) 

k=l i=l k=l j=l 
1 n 

-(J2(n-k + l) 2 + 2 (n-j + l)(n-i + l)(ei- ej )) 



k=l l<i<j<n 



□ 



Definition 6. As above, let us denote the average of the set {|| cp \\ 2 : P G V n } as 
|| c„ || 2 , the average center of mass ofV n . 

Lemma 5. || c„ \\ 2 = — - — . 

II n || ^ 

Proof. From Lemma [4], we have that for any P G? n , 
1 " 

II c P || 2 = —(Y^(n- k + 1) 2 + 2 ^ (ri-i + l)(n-i + l)(e i -e J -)). 

fc=l l<i<j'<n 
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We will replace (e* • ej) with r n to find || c n || 2 : 
1 " 

_(^( n _A; + l) 2 + 2 (n-j + l)(n-i + l)r n ). 

:j<n 
n j-1 



C II 2 

'-'77. 



fc=i i<*<i<™ 



i,£> - * + 1) 2 ) + 1 E E (<» - i + D(— < + 1) (^i)) 

/,• i j=2 i=i v ' ' ' 



Simplifying the first term, we have 

2n 2 + 3n + 1 



n 2 6n 

k=l 



Likewise the second term, 

n j-i 



2 

n 2 

j=2 i=l 

is simplified to 



^((n-j + l)(n-i + l) (wh)) 



1 n 

— — (f - ( 3 ™ + 4 )i 2 + O 2 + 7n + 5) j - (2n 2 + 4n + 2)) . 



n 2 (n — 1) ^ 



The fact that we are summing from j = 2 rather than j = 1 complicates this 
calculation. To fix this, we change the sum so that it starts at j — 1, and at the 
end, we will subtract off the j = 1 term: 

1 - 3n - 4 + 2n 2 + 7n + 5 - 2n 2 - An - 2 = 0. 

Because the j — 1 term is 0, we may change the summation to start at j = 1. 
This allows further simplification, and the second term is finally 

' 3n 4 + 2n 3 - 3n 2 - 2n 
12n 2 (n- 1) 

Combining these two terms, we have 

2 2n 2 + 3n+l 3n 4 + 2n 3 - 3n 2 - 2n 
" ° P " ~ 6^ 12n 2 (n - 1) 

(n + 1) 
12 



□ 
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3.5 Squared Radius of Gyration 

With the average center of mass defined, we may ask for the average squared 
difference between the center of mass and the vertices. This average difference is 
precisely the average squared radius of gyration. 

Definition 7. Let P G V n with vertices v±,...,v n . With cp as above, define the 
squared radius of gyration of P to be 



k=l 



Definition 8. We will define the average of the set {Rq(P) : P G V n } to be R 2 Gn , 
the average squared radius of gyration ofV n . 



Theorem 6. The average squared radius of gyration for all P G V n is R% n - 



n+1 

/few uj yyiuviun jui um c r n to -fi- G n — 

Proof. For any P G V n , we have that 



n 

k=l 

n 



v k f -2{v k -c P )+ || c P || 2 ) 

n z — ' 

k=l 

1 n n 1 

-T,\\v k || 2 -2(V^-c P ) + -(n || c P || 2 ) 

fc=i fe=i 
1 " 

- ^ II v k II 2 -2(c P • c P )+ || cp || 2 

fc=i 
1 n 

L E 



II 2 n n 2 

^fc - C P 

n '■ — ' " 

fc=i 

From Theorem |3 we have that the average value of || v k || 2 is fc< ^~^ ■ From Lemma 
[5] we have that the average value of || cp || 2 is ^1. Replacing these, we can find 
the average squared radius of gyration over all P G V n , that is, R@ n (rather than 
Rq(P) as above.) 



/ 1 k(n-k )\ n + 1 

K Cn - \ n 2^ n _l 



fc=l 

,2 



12 



n — 1 n + 1 n 1 
12(n- 1) 12 " 12 + 12' 



□ 
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Therefore the average squared radius of gyration scales like j|, agreeing with 
Zimm and Stockmayer's estimate |27j . 

3.6 Average Squared Radius of Gyration of Subsegments 
of Length k 

The scaling of the average squared radius of gyration of ideal rings and chains is 
influenced by the presence of knotting [16]. As a consequence, one expects that 
the average squared radius of gyration of subsegments will also be affected by the 
presence or absence of knotting. 

Determining a correlation between knotting and shape characteristics will be 
affected by the nature of knotting. If the knotted sections of the rings are small, 
involving relatively few edges, we may not see the influence of those edges in the 
squared radius of gyration, which can be obscured by the behavior of the unknotted 
complement. In order to detect local behavior, we study characteristics of shape 
that are local. 

In open chains, the average squared radius of gyration of subsegments of length 
k inside of larger chains of length n is identical to the average squared radius of 
gyration of segments of length k. That is, the ambient chain has no effect of the 
shape of the subsegment. This is not true for ideal rings, where we will see that the 
closure constraint plays a pivotal role in local shape. However, these characteristics 
may play a pivotal role in identifying and characterizing the knotted portions of 
open chains. 

Definition 9. Let P e V n . Define Pi± to be the translated subsegment of P of 
length k beginning with the i edge vector, e^. That is, Pi^ is a segment starting 



at the origin, where the j vertex is given by > e m . For ease of notation, we will 



So Pi } k is isomorphic to a subsegment of P, though we've done some relabeling. 
Now we would like to find the center of mass and the squared radius of gyration 




j 

relabel the position and edge vectors so that the j^ 1 vertex is given by v'j = e' r 

m=l 




we define the squared radius of gyration of P^ to be 



ii 



For P G V and k < n, define the average squared radius of gyration of subseg- 
ments of length k as 



1 



n 

3=1 

Now, we have the analogous definitions when we take the averages of cp i k and 
R 2 Gyk {P) over all P G V n . 

Definition 11. Let c k ^ n be average center of mass of a translated subsegment of 
length k, the average of {cp ik : i < n and P G V n }. For some n, let R Gnk be the 
average squared radius of gyration of subsegment of length k to be the average of 
{R Gik (P):PEV n }. 

As in the last section, in order to compute R Gnk , we must first find c k , n - 

„ , ll2 (2k 2 + 3A; + l)2n - 3k(k + l) 2 

Lemma 7. c k , n = — — . 

12k(n — 1) 

Proof. Let us begin with some || cp mk || 2 . 



3=1 3=1 
k j k j 



j=l i=l j=l i=l 

= ^&- i + 1 ) 2 + 2 ED fc -j + 1 )( fc - i+1 )( e '.-4) 



k k j—1 

J=l j=2 i=l 

Because {e[ -e'j) = (e' i+m ■ e^ +m ), the product is still independent of % and j, and 

we can replace (e- • e') with r n = — . This will let us take the average over all 

J n — 1 

such P mk to obtain || c kn || 2 : 



l -Aj2{k-r + lf + 2j2Y:{k-J + l){k-r + l) (^) 

\i=l 7=2 z=l V ' , 



i|2 
Cfc,n || 



There are two sums to evaluate. The first is straightforward: ^(/c — % + l) 2 = 

i=i 

: + 1)(2fc + ^ The second sum is 5 = 2 ^(jfe - j + l)(k - i + 1) (^y) • 

7=2 i=l \ n J 



12 



s = — Lj^fc-j + ixfc-i + i) 

n — 1 

j=2 i=l 
-1 * 



— £ (j 3 - (3* + 4)j 2 + (2k 2 + 7k + 5) j - 2(* + l) 2 ) 



n 

J=2 



Evaluating j 3 - (3A; + 4)j 2 + (2k 2 + 7k + 5)j - 2(fc + l) 2 at j = 1 we have 
1 — 3/c + 4 + 2/c 2 + 7/c + 5 — 2A; 2 — 4k — 2 = 0. Thus we can replace the lower bound 
of our sum, j = 2 with j = 1 with no penalty. 



-1 fc 

5 = V (f - (3k + 4)f + (2k 2 + 7k + 5)j - 2(k + l) 2 ) 

77,-1 *-~< 
3=1 

2k + 3k 2 - 2k 3 - 3k 4 
12(ra - 1) 

Combining these terms, we have 

2 _ 1 jfc(jfc + l)(2fc + l) 2k + 3k 2 - 2k 3 - 3k 4 
l|Cfc ' n|1 ~ P ( 6 + 12(n - 1) } 

(2k 2 + 3k + l)2n - 3k(k + l) 2 
12Jfe(n - 1) 

We note that when k = n, we have 

(2n 2 + 3n + l)2n - 3n(n + l) 2 



c II 1 



\2n(n - 1) 
n 3 — n n + 1 



12n(n-l) 12 



□ 



which agrees with Lemma [5] 

Theorem 8. TTie average squared radius of gyration of a subsegment of length k, 

_ (fc 2 -l)(2n-fc) 
',n,k ~ 12fc(»i-l) ' 



taken over all such subsegments in all P G V n is Rq , — — — 1 ^ 2 " k ^ 
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Proof. For some we have 

1 k 

i fc 



i=i 



i=i 
1 k 



4 



«i ii - ii c ^, fc 



From Lemma 



ll9 „ ll2 (2A; 2 + 3A; + l)2n - 3k(k + 1) 
c Pi fe |p is, on average, || c fcjn || = 



12Jfe(n - 1) 

Likewise, average value of || v'j || 2 , the end to end distance of the segment P^, is 
given in Theorem p| By replacing || cp ik || 2 and || v'j || 2 with the averages for these 
values, we can fina i?|. nfc : 



1 A j(n - j) (2k 2 + 3A; + l)2n - 3k(k + l) 2 
k^ n-1 Ylkin - 1) 

{k 2 -\)(2n-k) 
Y2k{n - 1) 



We note that when k = n, we have 



P 2 



□ 



(n 2 -l)(n) 
12n(n - 1) 
n + 1 
12 



which agrees with Theorem [8j 

Now for a fixed n, we have a function that returns the average squared radius 
gyration of a subsegment of length k. 
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3.7 Comparison Between Ideal Rings and Open Ideal Chains 

For comparison, consider ideal chains. From Lemma |4j which also holds for ideal 
chains, we know that for any W G W„, the squared center of mass is given by 



cw || 2 

n 



1 7 Jj2(n-k + l) 2 + 2 (n-j + lXn-i+lXei.e,)) 

\fc=l l<i<j<n / 



When the average is taken over all ideal chains W G W n , define the average cw„. 
Likewise, let the squared radius of gyration of the ideal chain W be R G (W), and 
let the average over all such ideal chains be RQ Wn . 

Lemma 9. Rcw n = 

Proof. First, we find || cw n || 2 , which is simplified by that fact that in W n , the 
average value of (e 4 • ej) is for i ^ j. 



~-i<j<n 



I (£>-*+i)'+2 y. 

\k=l l<i<j< 

K&- fc+i)2 ) 



n 

2n 2 + 3n + 1 
6n 

As in the first part of Theorem |8j the squared radius of gyration of some open 
chain, R G (W), is 



n 2 n n 2 

v k || - || C W || 

n A — ' 

k=l 



Replacing || cw || 2 with the average || c>v n || 2 and || t>& || 2 with its average value, k, 
from Lemma [TJ we can find the average radius of gyration for open chains, R GWn : 



1 

R G,w n = ~ ^2 k ~ I' Cyv " 

k=l 



1 /n(n + l)\ 2n 2 + 3n + 1 
n \ 2 / 6n 
n 2 -l 
6n 



□ 
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For open chains, the squared radius of gyration of a subsegment of length k is 




Length of Subsegment k 

Figure 3: Here we compare the average squared radius of gyration of subsegments 
of length k for open chains (blue) and ideal rings (red and dashed). 

As with squared end to end distance, we can see that the average squared 

radius of gyration of subsegments in an ideal ring is radically different from the 

average squared radius of gyration of subsegments in an open ideal chain. For 

ideal rings and chains of length n, in order to have the averages squared radius of 

gyration of a subsegment of length k for ideal rings to be within ^ of the squared 

radius of gyration of a subsegment of length k in an ideal chain, we must have 

each other, we must have that the length of the subsegment considered, k, is less 

than 5r . So for n = 1000, k must be 20 or less. 

Recall, with squared end to end distance, for the averages to be within 1% of 

th 

each other, the segments had to be less than ^ of the length, so k had to be 
length 10 or shorter if n = 1000. This suggests that the squared end to end distance 
characterization is more sensitive than squared radius of gyration of subsegments, 
as it is more affected by the closure condition. Once again, we observe that the 
local scale of an ideal ring differs significantly from that of a ideal chain at all but 
the smallest length scales. 
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4 Experimental Methods 



Generating a random sample from V n is much trickier than generating a random 
walk. With the latter, we can uniformly sample points edge vectors on the unit 
sphere. The complication in V n is the closure constraint: we cannot just generate 
random walks and hope that they close. 

In order to randomly sample V n we have two step process called the hedgehog 
method pQ. First, we find a starting point, then we use crankshaft rotations to 
sample the space. After some number of moves, our sample is independent of the 
starting point, and is a random element of V n . 

4.1 Hedgehog Method 

The hedgehog method begins with selecting n points uniformly on the unit sphere, 
and label them ex, ...,e n . Then, we add to that list each e^'s negative, — e^. Thus 
we have 2n edge vectors, e\, — ei, e^, — e 2 , e n , — e n . 




Figure 4: The edge vectors ei, — ei, e 2 , — e 2 , e n , —e n are plotted above, hence this 
is called the hedgehog method. 

We randomly permute these, getting a list of edge vectors e[, e' 2 , ■■■e 2n - Adding 
in each edge vector's negative ensures the all important closure condition is met. 

4.1.1 Crankshaft Rotations 

This is a good starting point, but for each edge vector, its exact opposite is also 
an edge vector, which is undesirable. 
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We finish the hedgehog method by performing crankshaft rotations. We begin 
a crankshaft rotation by randomly selecting two non-parallel edge vectors ej and 
e&. These are rotated about the axis determined by ej + e^, by a random angle 9. 
The form of these rotations is given by the following: 



6j + e& e,- e& . . x e^ . . . 

ej i — + — cos(#) + ^- j|Sin(0) 

i i || + e& || 

Because the sum ej + e^ is conserved, the modified sequence of edge vectors 
and vertices still satisfies the closure condition. 

Because any equilateral polygon can be deformed by a finite sequence of crankshaft 
rotations to the regular polygon, a finite series of crankshaft rotations will take 
us from one polygon to another random polygon pQ . We performed 6n crankshaft 
rotations on each sample polygon. 



5 Experimental Results 

We use the above methods to randomly generate ideal rings. Numerically, we will 
have two primary foci: how quickly does the average for an ensemble converge to 
these theoretical values, and how will knotting affect these characteristics? 



5.1 Convergence 

For various population sizes, 10 1 , 10 2 , 10 3 , 10 4 and 10 5 , we generated that many 
50 edged polygons with 150 crankshaft rotations. For each polygon, we found 
the squared radius of gyration, and took the average over the population. For 
each population size, we did this 10 times. We then compared these averages to 
the theoretical average. We define E(n) to be the difference between the exper- 
imental average and the theoretical average with n samples. Figure [5] shows the 
convergence result. 

For population sizes of 10000 and larger, we have excellent convergence to the 
theoretical average. Using linear regression, we can estimate that E{n) ~ CViT 1 ' 559 , 
providing for excellent convergence for large n. 



5.2 Knotting 

Knotted ideal rings, specifically trefoils, were generated by first constructing ideal 
rings of length 50, then calculating their knot type, and saving those of the given 
knot type. 
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Figure 5: For population sizes 10 1 , 10 2 , 10 3 , 10 4 and 10 5 we computed the average 
radius of gyration, and plotted it against the theoretical value for the same length. 
Note that this is a log scale plot of the data. 
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By the length of a knot in a ring we mean the length of the shortest subchain 
that contains the knot. The knot length is determined by examining subsegments 
of progressively longer length, starting at all possible locations. Each open segment 
is situated inside of a large sphere. A random collection of points on the surface of 
the sphere is selected. The ends of the segment are closed to each point, and the 
knot type of each of these closures is then calculated. This yields a spectrum of 
knot types, as in [UJ [T5J [T51 [TJJJ [23] • A segment is considered to be a trefoil if the 
closure is a trefoil with some tolerance (greater than 50%). The knotted portion 
is identified as the shortest segment which is a trefoil, for which the complement 
is unknotted. 




Figure 6: Here we compare the average squared end to end distance for a phantom 
population of ideal rings (blue) and a population of randomly sampled trefoils (red 
and dashed). 

Above we compare the average squared end to end distance for a phantom 
population of ideal rings and a population of randomly sampled trefoils. We can 
see that for length 50, the average for the trefoils is smaller than the average for the 
whole space. This suggests that for length 50, knotting compresses the polygon, 
making vertices closer together. 

Looking at the maximum squared end to end distance, we may ask what length 
curve has the same maximum, 10.2291. Solving, we have that a curve of 39.89, 
approximately 40, has the given maximum. That suggests that the average short- 
ening caused by knotting is about 10 edges of length, and that on average a trefoil 
of length 50 has an effective length of 40. 
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The average length of the knotted portion of these trefoils is 16.4. We predict 
that the difference between the length 10 reduction we saw above and the average 
trefoil length, 16.4, can be accounted for by examining the end to end distance 
of the knotted portion. As in Figure [7], the difference could be explained by an 
average knotted portion of 16 to 17 edges with an average end to end distance 
around 6. 




Figure 7: A hypothetical average knotted section with knot length 16 edges and 
end to end distance, (red and dashed), of 6. 




Length of Subsegment k 

Figure 8: Here we compare the average squared radius of gyration of subsegments 
of a phantom population of ideal rings (blue) and of a population of randomly 
sampled trefoils (red and dashed). 
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As with squared end to end distance, the average for the trefoils is smaller 
than the average for the whole space. Again, this suggests that for length 50, 
knotting compresses the polygon, making vertices closer together. We know that 
ideal trefoils, and indeed any collection of some fixed knot type, will for some small 
number of edges have mean squared radius of gyration less than the the mean for 
the total population, and for larger number of edges, they will have a greater 
average squared radius of gyration J5J. For example, the average squared radius 
of gyration of trefoils is smaller than the average squared radius of gyration of the 
whole population for lengths less than 175, and the opposite is true for lengths 
greater than 200 [5]. We expect similar behavior for squared end to end distance. 

We can set the average squared radius of gyration of the population of trefoils, 
3.5768 equal to our function for squared radius of gyration, Solving for n 

we have n = 41.9216, suggesting that the average shortening caused by knotting 
is about 8 edges of length. We can compare this with the average shortening 
prediction from squared end to end distance, 10. These differ by 4% of the total 
length, 50. 

6 Conclusion 

These theoretical averages have many potential applications. Primarily, they can 
be used as a criterion to determine effectiveness of sampling methods. By com- 
paring the the squared radius of gyration of subsegments or the squared internal 
end to end distance of polygons generated a given sampling method and the the- 
oretical average, we can determine the effectiveness of the generation scheme and 
corroborate numerical simulations. Further, in the above generation method, we 
can use this convergence to determine how many crankshaft rotations are needed 
to sample the space of polygons uniformly, as in pQ. 

As the previous section highlights, squared end to end distance and squared 
radius of gyration of subsegments may be used to predict knot length, which is 
computationally expensive to calculate. This will allow us to examine the growth 
of the knot length as n — > oo, to determine if average knot length is bounded, or 
grows proportionally with n 1 / 2 or n, allowing us to ascertain if knotting is strongly 
local, local or global on average pfl [El US EEl ED] • 
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